Symmetrical charge-charge interactions in ionic solutions: 
implications for biological interactions 
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As is well known in electrolyte theory, electrostatic fields are attenuated by the presence of mobile 
charges in the solution. This seems to limit the possibility of an electrostatic repulsion model of 
biological interactions such as cell division. Here, a system of two charges in an ionic solution is 
considered. It is found that in the context of the symmetries of the system, the electrostatic repulsion 
between the two is considerably increased as compared to the electrostatic repulsion between two 
bare charges in a dielectric. This increase in repulsion, directly resulting from interactions between 
the symmetrical parts of the system, was found to be dependent on the magnitude of the charges 
and the separation between them. It was also found that this increases reaches a steady state 
for separation greater than a solvent determined length scale related to the Debye length. These 
findings strongly suggest that electrostatic interactions can play a crucial part in the physical forces 
that are involved in biological interactions. 

PACS numbers: 87.17.Ee, 87.10.Ca, 87.17.Aa, 87.50.cf 
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INTRODUCTION 

Biological cell division is the corner-stone of life. From 
a practical stand point, its malfunction is responsible for 
many diseases. Cancer and neural regeneration are two 
fields of research relevant to countless people that can 
gain tremendous advancements from a basic understand- 
ing of the physical mechanism of cell division and an abil- 
ity to manipulate them. Even viral and bacterial infec- 
tions require the division of cells and the propagation of 
genetic information. Because of cell division's primordial 
necessity for life, the fundamental underlying processes 
of it must be general. Because of these fundamental val- 
ues, it is a candidate for systems where general biological 
interactions may be observed. 

All life shares DNA as the basic code of its partic- 
ular form. All life also shares the twenty amino-acids 
that make up proteins, the matter of life. This points 
in the direction of a common ancestor. It is also plausi- 
ble to assume that the mechanisms for DNA segregation 
will rest on the same fundamental physics across differ- 
ent life forms. Probably the first thing the first living 
cells did was divide. Hence, it is required that the phys- 
ical processes were simple, and did not involve complex 
biological machinery that would only develop later. If 
life originated by the chemicals of some ancient sea, act- 
ing freely under the forces of nature, those primordial 
cells were some form of dividing DNA. An extraordinary 
process. 

There are four known fundamental forces. Two of 
them, the strong and weak nuclear forces, have to do 
with interactions on the atom's nucleus length scale. Per 



current knowledge, these forces do not contribute to life 
beyond their role in the creation of matter. A third force, 
gravity, has to do with large length scales, it is negligi- 
ble in the biochemistry of living organism. This leaves 
us with the electromagnetic force as the sole interaction 
responsible for the differences between living and non- 
living matter. In fact, quantum electromagnetism is the 
basis of atomic structure and interaction. When moving 
from a few atoms to the many atom systems prevalent in 
living organisms it is found that quantum effects become 
mostly negligible, and classical electromagnetism (EM) 
is used. Many models use EM to describe biological pro- 
cesses. In the realm of cell division, some form of EM 
must hide in the complexity. 

While it is known that tubulin, a protein related to 
cell division, posses a permanent dipole moment [3- 3] , 
and that electric charges play a role in cell division 
H, 12 1, other types of models were also proposed. [fj|-[ll| 
A major objection to electrostatic models of cell division 
is Debye screening for ionic solutions. 0, [l3[ Considering 
the fields of a single charge in a solution of freely mov- 
ing charges, one finds that the electrostatic potential is 
attenuated by an exponentially decaying envelope due to 
the presence of the mobile charges. The assumption be- 
ing that the mobile charges will be distributed according 
to the Boltzmann distribution. Since the order of magni- 
tude for the Debye screening length for biological matter 
is 1 nm, it was not expected that electrostatic repulsion 
can play a role in cell division. 

Here it will be shown that analysis of a symmetrical 
system, as is the case for example for the two kinetochores 
attached to the genetic code, shows an opposite trend. 
While the individual electromagnetic interaction is expo- 
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ncntially attenuated as prescribed by Debye screening, 
the many-body force of repulsion between the two sym- 
metrical parts is increased as compared to bare charges. 
It is the intent of this work to argue that EM models can 
be constructed to describe cell division and other biolog- 
ical interaction. As briefly discussed above, there is no 
other candidate for a fundamental force except EM. 



MODEL 

Consider the Poisson equation for the potential $, in a 
system of a point charge, qi , at the origin, with stationary 
charged background (charge — q m ) and a cloud of mobile 
charged particles (charge q m ) 

V 2 $ = -47TQ 1 5(x) - 4^ m n [e-^*/ fcBT - 1]. (1) 

x is the vector position, tiq the particle density, and ksT 
the thermal energy. For a background of this equation 
see Jackson. 13[ The second term on the right hand side 
of Eq. |T]) comes from the Boltzmann distributed mobile 
charges and the stationary ones. 

It is instructive to first examine charge conservation 
and the gauge invariance of the electrostatic potential. 
We need to separate two cases: a finite or infinite sys- 
tem. For a finite system a partition function is needed 
to normalize the exponential term in the Boltzmann dis- 
tribution. This partition function will guarantee both 
the gauge invariance and conservation of charge. In the 
limit of infinite distances, the perturbation of the back- 
ground/mobile charges due to the point charge should 
become negligible and hence 
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This leads to the standard boundary conditions ^x^oo = 
0. A partition function normalization is possible if one 
wishes to add a constant to $. 

If one assumes that the EM energy is small compared 
to the thermal energy, one can approximate the term in 
the square brackets in Eq. (JTJ as linear in $ and the 
resulting differential equation leads to an exponentially 
decaying Coulomb field with a length scale of 1/e decay 
given by the Debye length. Since the Debye length is 
of the order of a few protein helix turns it was judged 
that electrostatic repulsion is of negligible importance in 
biological cell division. 

The crucial point that this work picks up on is that the 
dividing biological cell is inherently a symmetric process. 
This implies that we are interested in the interaction be- 
tween two charges in such an environment. As will be 
discussed in a future, more biologically oriented publi- 
cation, the existence of the two charges can come about 



from the dipolar nature of tubulin. Here the focus is on 
the physical interaction of a system as depicted in Fig. [T] 
A mirror symmetry is taken between the two parts of the 
system. 





FIG. 1: Schematic representation of the model. Two charges, 
qi and q%, are separated by the colored plane in the figure 
which is perpendicular to the line joining the charges and 
placed equidistant between them. This plane defines two vol- 
umes, vi and V2 as shown in the figure. 

The modified Poisson equation for the two charges, 
9i = qi = q immersed in an aqueous solution with dielec- 
tric D, is given by 



D-V 2 $ = -47rg[<5(x+x )+(5(x-xo)]-47rg„ 



n e 



-q m <S>/k B T_ 

(2) 

taking a coordinate system where the charges are at ±xo 
as described in Fig. [TJ As before, if \q m $/Dk B T\ « 1 
one can approximately solve Eq. © and obtain 



q g-fcoUx— x || g-fc^Hx+xoll 

= — \—r, n- + -n n~l 



(3) 



with kjj = Annoq^ / DksT the Debye constant in a di- 
electric. 

This is the field between the two charges, attenuated 
by Debye screening. For a homogeneous system the elec- 
tromagnetic fields far away from the two charges is at- 
tenuated. Here the interest is in the force of separation 
between the charge densities of the two volumes delin- 
eated by the symmetry of this problem. As displayed in 
Fig. [T] these volumes are separated by the plane contain- 
ing the origin and perpendicular to the line joining qi 
and 172- The Coulomb force between the volumes V\ and 
Vi can be formally written as 



F s = 



P(x)p(xQ 
xeVi Jx'ev 2 D ||x — x'|| 2 



d'V d 3 x, 



(4) 
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with 



p(x) = q[S(x + x ) + 6(x - x )] + q m n [e 



-q m <S>/k B T 



The integral in Eq. Q can be re-represented to aid 
in its interpretation. Integrating out the delta functions, 
defining 7(x) = e -9m*(x)/fc fl T _ ^ anc j usm g ^ ne symme- 
try of the system F s becomes 



F h + F„ 



F r 



with 



F r 



F a 



F b = 



2n Q q m q 
D 



D||2x || 2 



7 (x) 



■d^x 



x 



(5) 



(6) 



(7) 



D 



xeVi Jx'ev 2 Il x - x 



7 (X)7(X 2 d 3 xdV (8) 



Equation (JS|) breaks up the interaction forces into three 
components. Ff, is the Coulomb force between the two 
charges external to the ionic solution. F a represents the 
force between the induced ionic charges and the bare 
charge in the opposite volume. Specifically, it represents 
the interaction of charge 1 with the ionic cloud surround- 
ing charge 2 plus the interaction of charge 2 with the 
ionic cloud surrounding charge 1. Finally, F r represents 
the force between the induced ionic charges in the two 
volumes. Since EM is a linear theory, the total force 
between the two volumes is the superposition of the in- 
dividual forces. The symmetry of the system considered 
dictates that both Ff, and F r are repulsive, while F a is 
attractive. 

Up to now the derivation for the two particle system 
is general in terms of the relationship between thermal 
and EM energies. In principle, if one can solve Eq. @ 
and obtain 7, one can use Eq. to obtain the force be- 
tween the charge distribution in the daughter cells. One 
should note that the general solution is not symmetric 
with respect to charge. This comes about since for an 
infinite system we are unlimited in the amount of mo- 
bile charges we can add. Positive charge density, on the 
other hand, is introduced by the removal of this mobile 
charges. Hence, it is bounded by the initial, unchanging 
density of the immobile ones. 

In practice, however, it is well known that biologically 
relevant ionic energies are significantly smaller than ther- 
mal ones. One major reason for this is the abundance 
of water molecules. Due to their polarity these tend to 
aggregate around charges and limit the access distance 



for ions. If one assumes that the ionic energy is small 
compared to the thermal energy, one can approximate 
7 as —q m $/k B T. For convenience, define from Eq. ([3]) 
%jj = DQ/q, this enables us to write for F a and F r : 



F a = -2 



n {q m q) 2 
D 2 k B T 



xo 



F r 



D 3 (k B ry 



(9) 



(10) 



xGVi Jx'ev 2 



In this regime, changes in the charge density are small, 
and as exhibited by Eqs. (|9|10|) are symmetric with re- 
spect to charge reversal. If the charge density fluctua- 
tions are larger this symmetry is lost, as mentioned pre- 
viously. We should now consider the finite size screening 
of the EM forces near the bare charges. Since this in- 
volves the collective behavior of a large number of water 
molecules it is an entropic effect and will depend on the 
temperature. For convenience it will be assumed that 
q m $ < 0.25 *k B T, with T = 300K. This condition is also 
necessary to satisfy the requirements of the linear approx- 
imation for 7. A physical interpretation for this condition 
is a layer of water of radius a = qq m /2 ■ 0.25Dk B T sur- 
rounding the charges. This layer inhibits closer approach 
between the various charges. Thus, to the edge of the 
water layer, EM fields propagate through an ion-less di- 
electric. That is, Eq. Q should be modified to separate 
distances smaller and larger than a. Specifically, the EM 
field decays exponentially only for distances greater than 
a and a continuous field is assumed at a. This amounts 
to subtracting a from the exponent in Eq. (|3]) for the 
potential outside the dielectric region. 



RESULTS 

The analysis of the treatment outlined above was done 
using trapezoidal numerical integration coded in FOR- 
TRAN. Given on the y-axis of Fig.[2]is the total force be- 
tween the charge distributions in the two daughter cells, 
F s , normalized by the Coulomb force between two core 
charges taken as protons in a water dielectric, Fb. On 
the x-axis of Fig. [5] the distance between the charges is 
given. We see that up to a separation of approximately 
3 nm F s « Fb- This is due to the limit we put on the 
distance between the charges which for this case is given 
by 2.7 nm. As the separation is increased, effects from 
the ionic field become more noticeable with ions in the 
two daughter cells repealing each other. This regime of 
sharp transition, where F s becomes almost an order of 
magnitude larger than Fh, reaches a plateau around 4 
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nm after which F s — » 6 • i<&. That is, due to the symme- 
try between the charge distribution in the two daughter 
cells the repulsive force between q\ and qi is increased 
six-fold as compared to their Coulomb repulsion and this 
ratio remains approximately constant over the range of 
the calculation. 
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FIG. 3: Ratio between F a and Fb, as a function of the sep- 
aration between the charges. In this case three q-values are 
presented. The previous case of q — le is barely discernible 
on this scale. The case of q = 6e reaches saturation ratio of 
about 300. The case of q — lie does not reach saturation 
since for this case a = 30 nm. 



DISCUSSION 



FIG. 2: Ratio between volume separation force, F s , and the 
force between two proton charges in a dielectric, Fb, as a 
function of the separation between the charges. We see that 
for ||xo|| < a = 2.7 nm the ratio increases as a the separation 
increases and it reaches saturation of around 6.5 for ||xn|| ~ 
4 nm, corresponding with a. The curve through the data 
points is used merely to aid the eye. 



Results for q — 1, 6, 11 e with e the charge of the proton 
are given in Fig. [3J Again, we see that for small sepa- 
ration F s w Fb. As the separation is increased, contri- 
bution to the repulsive force from the ionic field become 
more notable. These are most likely due to contributions 
from the cleavages in the misaligned spheres defined by 
a around q\ and q2- For q = lie, a = 30 nm, explain- 
ing the inability of this case to saturate by | |xo | =20 nm, 
the largest calculated separation. For the case of q = 6 e, 
a = 16 nm, which corresponds well with the saturation 
point for that curve. Note that the boundary condition 
for the numerical integration were set at 31 nm due to 
computational reasons. The proximity of the boundaries 
causes the decay observed in Fig. [3] for the q = 6e case 
past saturation. A similar effect was observed for q = le. 
Note that the single point in Fig. [2] at ||xo|| = 30 nm 
was calculated with the outer boundary at 61 nm. This 
boundary effect also critically limits the accuracy of the 
results for q — lie, and is exhibited by the crossing of 
the lines for q = lie and the q — 6e at a separation of 
approximately 8 nm, due to the greater a of the q = lie 
case. 



A system of two charges in an ionic solution was con- 
sidered. It was found that in the context of Debye ionic 
screening the two symmetrical parts delineated by these 
charges experience an effect that increases the electro- 
static repulsion between them as compared to the ion- 
free case. This increase in repulsion is dependent on 
the magnitude of the charges and the separation between 
them. It was also found that this effect reaches a steady 
state for separations greater than a solvent determined 
length scale related to the Debye constant. These find- 
ings strongly suggest that electrostatic interactions can 
play a crucial role in the physical forces that are respon- 
sible for biological interactions such as cell division. It is 
also interesting that other deviations from electrostatic 
interactions were observed, for example for the system of 
lipid membranes. 
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A separate question, under current investigation, is 
that of the mechanisms and processes that would be 
able to reproduce the observed symmetrical configura- 
tion of the dividing biological cell, and their relationship 
to the configurations discusses here. The main ingredi- 
ents of this formulation are: dipolar tubulin, symmetry in 
its connections to the centrosome, and the kinetochore. 
Since the tubulin connects with a preferential polarity 
to the centrosome, through the dipole-dipole interaction 
along the tubulin fiber, charges of identical polarity will 
connect to each of the kinetochores. These two repuls- 
ing charges in the kinetochores could be the basis of the 
separation of DNA and biological matter. 
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